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RESPONSE OF AN INFINITE ELASTIC PLATE TO 
AXISYMMETRIC INITIAL VELOCITY DISTRIBUTIONS WITH 
APPLICATION TO HYPERVELOCITY IMPACT* 

By Robert J. Hayduk 
Langley Research Center 

SUMMARY 

Classical plate theory is used to study the response of an infinite plate to an axi- 
symmetric initial velocity distribution over a central region. The expression for the 
bending stress is reduced to a single definite integral which depends upon the initial veloc- 
ity distribution. The deflection is reduced to the Hankel inversion integral for three 
example distributions — conical, cylindrical, and spherical. The analysis indicates that 
the conical distribution causes larger stresses at the origin and deflections of the plate 
at short times after initiation than either the spherical or cylindrical distribution. 

A simple relationship is derived for the failure thickness of an elastic plate 
responding to an axisymmetric initial momentum, monotonically decreasing in the radial 
direction. The failure -thickness relationship is shown to correlate with available exper- 
imental hypervelocity impact data where the initial momentum was induced by disk pro- 
jectiles of various diameters. 


INTRODUCTION 

This paper presents a study of the response of an infinite, thin, elastic plate to an 
axisymmetric initial velocity distribution over a central region. The problem is directly 
related to an impulsively loaded plate, where the impulse is sufficiently distributed to 
ensure that the elastic plate action is the primary response and that the boundaries are 


*Part of the information presented in this paper was included in a thesis entitled 
"The Response of a Single Wall Space Structure to Impact by Cometary Meteoroids of 
Various Shapes" which was submitted in partial fulfillment of the requirements for the 
degree of Master of Science in Engineering Mechanics, Virginia Polytechnic Institute, 
Blacksburg, Virginia, June 1968. 



so remote that they do not influence the response until some time after the maximum 
stress occurs. 

Classical plate theory, which neglects transverse shear and rotary inertia, has 
been used by many investigators to study the response of plates to impulsive loads. Ref- 
erence 1 contains solutions to impulsively loaded infinite plate problems frequently ref- 
erenced in the literature. Particular problems are an impulsive point force at the origin, 
an impulsive force uniformly distributed over a circular area, and an impulsive Gaussian 
distribution of pressure. In reference 2 the interaction of the spray from a meteoroid- 
penetrated bumper with the main wall of a double-walled structure was analyzed as a 
Gaussian distribution of initial momentum impacting on an infinite plate. Reference 3 
showed that the deflections at a distance from the impact point of a large plate centrally 
impacted by a bullet agreed after a short time with those obtained by mathematically sim- 
ulating the impact by an impulsive point force of finite duration on an infinite plate. Also 
presented in reference 3 are solutions to the impulsive point force problem for various 
time functions: finite duration, step function, and unit impulse function. 

A large, highly concentrated impulse on a plate will cause plastic shear deforma- 
tions. If, however, the impulse is distributed over a sufficiently large area, bending 
stresses will prevail, and the response of a plate can be well estimated by classical elas- 
tic plate theory. Axisymmetric distributed impulse loadings occur under certain circum- 
stances on the second wall of a bumper -protected structure and on a single -wall structure 
impacted by cometary meteoroids of sufficiently low density. In the first case, the 
momentum of the spray from the pulverized bumper and meteoroid material imparts an 
initial velocity to the main wall; in the second case, the momentum of the low -density 
low-strength cometary meteoroid imparts an initial velocity to the single wall. 

A general formula is derived for the maximum plate stress, and the response to 
three particular initial velocity distributions is examined in detail. A plate failure - 
thickness relationship is derived for an axisymmetric monotonically decreasing momen- 
tum impulse. Cometary meteoroid impact on a single-wall structure is considered as a 
momentum impulse and the form of the failure-thickness relationship is compared with 
available hypervelocity impact data. 


SYMBOLS 

a radius of central region given initial velocity 

A(p), B(p) functions of p 
B(o-,i3) 


2 


beta function where a and 8 are the arguments 

r(os + (3) 




d 

E 

f 

1 F 1 (a;^;4), 

h 


diameter 

modulus of elasticity 
geometric constant 

F 2 (at ; /3 ;y;£) hypergeometric functions where 

arguments 

plate thickness 


and 


i-fl 

Jg(p? 7 ) zeroeth-order Bessel function of first kind 
J ^(P 7 ?) first-order Bessel function of first kind 
J 3/2 (?) fractional -order Bessel function of first kind 

k = . m 2) 

l length of projectile parallel to direction of motion 

m mass 

M magnitude of momentum per unit area at 77 = 0 

n summation index 

p Hankel transform parameter 

Re( ) real part of complex variable 

r radial coordinate 

s Mellin transform parameter 


are the 
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shape factor 


nondimensional variables of integration 

initial axisymmetric velocity distribution over a central region of radius a 

dimensionless function specifying distribution of initial velocity or initial 
momentum 

magnitude of initial velocity distribution (value at 77 = 0) or projectile 
velocity with subscript 


w( 7 ?,t) 


transverse displacement of midplane of plate 


distance measured from and normal to midplane of plate 


r(o-), r(j3) gamma functions 


biharmonic operator 


nondimensional radial coordinate, r/a 


V^ofov) kernel of Hankel transform 


Poisson's ratio 


mass density (plate mass density without subscript) 


stress in general or stress at origin of plate 


nondimensional time, t/k 
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Subscripts: 


CO 

conical 

cr 

critical 

cyl 

cylindrical 

f 

failure thickness 

max 

maximum 

P 

projectile 

sph 

spherical 

96 

circumferential 

9V 

radial 

V 

order 


Superscripts: 

\,v constant exponents 

A bar over a symbol denotes the Hankel transform; a tilde over a symbol denotes 
a nondimensional value. Dots over symbols denote differentiation with respect to time. 

GOVERNING EQUATIONS 

Definition of the Problem 

Consider a thin, elastic plate infinite in extent with coordinates and dimensions as 
shown in figure 1. The plate's response to an initial axisymmetric velocity distribution 
v(?]) over a central region of radius a can be predicted by classical plate bending 
theory. The initial displacement of the plate is zero and the boundary and continuity con- 
ditions are prescribed by the transformed governing equation, that is, zero slope, finite 
curvature, and finite shear per unit length at the origin, and zero deflection, slope, curva- 
ture, and shear per unit length at infinity. 
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Transformed Governing Equation 


The governing equation for the bending of a thin plate from linear, elastic, small- 
deflection theory is 

V 4 w = -k 2 (la) 

at 2 

where 

k = - M 2 ) (lb) 


and the operator is 


V 4 = 


1 d_ 

7 ] dr) 



1 a_ 

77 dr) 



(lc) 


The deflection w(r7,t) is the transverse displacement of the midplane of the plate, a is 

the radius of the central region which is given the initial velocity distribution, c = 

h is the plate thickness, p is Poisson's ratio, and 77 is the nondimensional radial 
coordinate r/a. 

Since the plate is infinite in extent, the problem is amenable to a zeroeth-order 
Hankel transform in the nondimensional radial coordinate 77. The transform and its 
inverse are, respectively, 



f(p) = j r/f (77) J 0 (pt?) d?7 

and 


(2a) 


fO?) = j Pf(p) J 0 (P*?) d P ( 2 b) 

where p is the Hankel transform parameter and 7 jJq(pt 7) is the kernel of the transform. 

Application of the Hankel transform to the governing equation yields the following 
intermediate equation after integration by parts several times 



1 _ 8 _ 

7? 817 



+ p if) J l( p77 * " p2r? 5 ^- j o(p^ “ P 3 ^!^)) + k 2 ^f-+ p 4 w = 0 

J 77=0 


( 3 ) 
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The term in braces vanishes at the upper and lower limit on 77 upon application 
of the appropriate boundary conditions at infinity and continuity conditions at the origin. 
At infinity, the deflection, slope, curvature, and shear per unit length are selected as 
zero and of the order necessary for the terms in braces to vanish. At the origin the con- 
tinuity conditions are selected as zero slope, finite curvature, and finite shear per unit 
length. If the time is nondimensionalized to T = p the transformed governing equation 
is 

^+P 4 w = 0 (4) 

dr 2 


GENERAL SOLUTION 


The general solution to equation (4) is 

w(p,t) = A(p) sin(p 2 -r) + B(p) cos(p 2 r) (5) 

where A(p) and B(p) are functions of the initial conditions on the plate. The function 
B(p) is zero since the initial displacement, and hence the transformed initial displace- 
ment, is zero. The first derivative of equation (5) with respect to t relates A(p) to 
the transformed initial velocity distribution; that is, 


^(p,t) = IdwJkri 

k dr 

w(p,t) = — p 2 A(p) cos(p 2 r) 
k 


(6a) 

(6b) 


Hence 


A(p) = kp -2 w(p,0) 

(7a) 

A(p) = kp" 2 v(p) 

(7b) 


Substitution of equation (7b) into equation (5) with B(p) = 0 and application of the Hankel 
inversion integral (eq. (2b)) yields the displacement of the plate 

P 00 

w (??, t) = k \ p _1 v(p) sin(p 2 T)J 0 (p77) dp (8) 
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Once the initial velocity distribution of particular interest is transformed, the displace- 
ment reduces to the inversion integral as shown by equation (8). The infinite integral 
usually defies analytical evaluation and must be evaluated numerically. Approximate 
methods for evaluating integral transforms and infinite integrals are available in standard 
mathematics references and usually exist as subroutines in computing programs. 


With the displacement in inversion integral form, the stress, strain, shear, and 
moment distributions throughout the plate are easily obtained by differentiation and appli- 
cation of the same approximate methods for evaluating the integrals. Only loadings mono- 
tonically decreasing in the radial direction are considered in this paper. For such 
loadings the maximum stress, strain, and bending moment occur at the center. For these 
maximum responses the solution can be reduced to a single finite integral which is a func- 
tion of the initial velocity distribution only. 


Since 


f=- =o, 

97 ?J?7=0 


lim I2w B 8 2 wl 


77 — 0 V d V drp 


-1?7=0 


( 9 ) 


and the functional forms of the stress, strain, and bending moment at the origin are iden- 
tical; therefore, the discussion in this paper is limited to the stress. From plate theory 
the normal radial and circumferential stresses at point z measured normal to the mid- 
plane are 


_ zE ( 9^w y. 9w\ 

^ a 2 (l - fj. 2)la??2 V dr] J 


and 


a ee ~ 


zE / 1 9 w , ,, 9 ^w 


— + u 

'vi r\ 1 


a^(l - /j,2)y? dV 0T7 2 
which, because of equation ( 9 ), reduce to 


= -zE afw 

r]=0 U ^r]=0 a 2 (l - n) d V 2 J v=0 




(10a) 


(10b) 


(ID 


Substitution of the transform integral (eq. ( 2 a)) for v(p) into equation (8) yields 


w(t?,t) = k f uv(u) C p -1 J 0 (pu) J 0 (p77) sin(p 2 r) dp du 
Jq Jq 


( 12 ) 
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if the order of integration is reversed. The upper limit on the first integral is unity 
because v( 7 j) = 0 for 1 < 77 < The second derivative of w(t],t) with respect to r\ 
evaluated at r\ = 0 is 


3 2 W 

9t7 2 


7J=0 


= - uv(u) y p J q ( pu) sin(p2r) dp du 


The infinite integral is evaluated in reference 1 as cos^i-j provided r ^ 0. The 
maximum stress occurs at the origin (eq. (11) with z = h/2) and is 


*< t ) = (5 ^ e p) 


1/2 


fX uv<u> cos (^) du _ 


(t / 0) (13) 


Note that the maximum stress is a function of time. Thus the maximum stress (strain 
and bending moment resultant per unit length) at the origin has been reduced to a single 
finite integral which can be evaluated for any particular velocity distribution of interest. 


SOLUTIONS FOR SPECIFIC VELOCITY DISTRIBUTIONS 


In the following sections cylindrical, conical, and spherical initial velocity distri- 
butions are considered to demonstrate the application of equations (8) and (13). In the 
applications section the arbitrary magnitudes (V terms) of the example initial velocity 
distributions are related to a general axisymmetric momentum impulse monotonically 
decreasing in the radial direction. 


Cylindrical Distribution 

Consider the case where the central portion of the plate is given a uniform initial 
velocity; that is, 


V C ylM = V cy i 


(0 < 7] < 

<! 
0 
<< 
t— * 

11 

0 

(1 < V < 

Vcyifa) transforms as 

cyl(P) = V cy i J VJ 0 (PV) d h 




(14) 


(15) 


V C yl(P) = V cyl p-^ifo) 
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( 16 ) 


Substitution of v cyl (p) into equation (8) gives the displacement of the plate as 

C°° 

w cy l(^ T ) = kV cyl J 0 P“ 2j i(P) J o(P^ sin(p 2 r) dp 
The maximum stress at the origin (eq. (13)) becomes 

,1/2 V 


a cy i(T) = Ep ) Iq 1 u cos {£) du (T ^ 0) (17a) 


which readily integrates to 


7 cyl (T) = |(^T^ L) Ep ) 1/2 V cyl sin ( 4T ) _1 (t ¥ 0) (17b) 


Conical Distribution 

This distribution is expressed as a velocity linearly decreasing with radius; that is, 


v co fa) = v co(! - V) 

Vco^) = 0 


(0 < 77 < 1) 

(1 < 77 < °°) 


Application of the transform to this v(? 7 ) results in two integrals: 
v C0 <P> = V co f v dri - V co J n ^ J 0 ^) 


(18) 


(19a) 


The first integral is identical to equation (15); the second is evaluated in the appendix. 
For the conical distribution 


v co (P) = Veop-ijjfp) - i V co 1 F 2 (!; 1 ,!;-ip 2 ) (19b) 

The displacement becomes 

c °° if 00 

w co (t 7 ,t) = kV C0 J p" 2 J 1 (p) J 0 (p7?) sin(p 2 r) dp - - kV c0 J p _1 X F 2 Jq(p??) sin(p 2 r) dp 

( 20 ) 

where the arguments of the hypergeometric function remain unchanged. Comparison of 
equations (20) and (16) indicates that the first integral represents the displacement due to 
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a uniform velocity distribution and the second integral represents the correction neces- 
sary to account for the linear decrease of velocity with radius. 

For this case the maximum stress at the origin (eq. (13)) becomes 

CTco(t) = Ep ) 7 .C u(1 " u)cos (i^) du (r ^ 0) (21a) 

Integration by parts several times and a change of variables yields the Fresnel sine 
integral 

a co( T ) = + ^ E Pj v co rl/2 ^ sin u2 du ( 21b ) 

For very small r, approximate expressions exist for the integral to facilitate numerical 
integration. 


Spherical Distribution 

The distribution is given as 


Vsphfr) = V sph (l - t? 2 ) 1 / 2 (0 < rj < 1) 

v sph(^) = 0 (1 < V < °°) 

The transformed v g pj 1 (?7) is 

^sph^ = ^sph " ^ 2 ) ^ J 0 (P^ 

Evaluation of this integral (see appendix) yields 

1/2 

^sphfc) = p“ 2 /2 Jg/gfo) 

The displacement for this example becomes 

w sph(’b T ) = kV sph(|) ^ P" 5//2 J 3/2 ^ J o(P^ sin(p 2 T) dp 


( 22 ) 


(23a) 


(23b) 


(24) 
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and the maximum stress at the origin 

CT sph( T ) = ( ^ff^ Ep ) 7 "^T J 0 U(1 " U2)1/2 C ° S (^) dU (T ^ 0) (25a) 

The integral is evaluated in the appendix and yields 


a M » 1/2 /3(1 + m) e A 1/2 „ V (-l) n (4T)- (21 

»sphW - — l-y-rf- E f] v sph 2, /■ — ? 

' ' n=0 r 2n + - 


-(2n+l) 


(t ^ 0) (25b) 


For very small r, an asymptotic expansion for the series facilitates convergence. 


RESULTS AND DISCUSSION 


As a basis for comparing the effects of the three example distributions, the velocity 
magnitudes (V terms) were selected to give the plates the same total initial momentum. 
This basis was chosen to facilitate application to a general axisymmetric initial momen- 
tum impulse monotonically decreasing in the radial direction. 

If the diameters of the central regions and the total initial momentums are equal 
in all three cases, the velocity magnitudes are related. The total initial momentums 
for the three cases are: 


Momentum^ = 7ra^hpV C yi (26a) 

Momentum co = ^ 7ra2hpV C0 (26b) 

O 

9 o 

Momentunigph = Tra^hpVgph (26c) 

The second and third momentum expressions equal the first if 

Vco = 3v cyl <2M > 

V sph = | V cyl < 26e > 


To demonstrate the effect that the shape of the initial velocity distribution has on 
the displacement of the plate, the deflection at t = 0.15, 0.4, 0.6, and 1.0 (eqs. (16), (20), 
and (24)) is presented in figure 2 for the equivalent magnitudes of V C yp 3V C yj, and 

V C yi for the cylindrical, conical, and spherical distributions, respectively. The 
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ordinate in the plots is the nondimensional displacement obtained by dividing the actual 
displacement by kV C yj. 

Figure 2 indicates that the conical initial velocity distribution causes the largest 
displacement of the plate at a particular time and the cylindrical distribution, the least. 

The differences in deflection caused by the various velocity distributions are large initi- 
ally, and then diminish with time. By t = 1.0, the deflection of the plate is nearly the 
same for all three distributions. 

The maximum bending stress at the origin of the plate for the three example velocity 
distributions (eqs. (17b), (21b), and (25b)) with the magnitudes given by equation (26) is 
presented in figure 3 for t ranging from 0 to 0.7. The ordinate parameter is the non- 

dimensional stress | * ~ LL — - 

\3(1 + p)Ep 

then the stresses rapidly converge. The peak stresses for the conical and spherical dis- 
tributions are 1.58 and 1.22 times that of the cylindrical distribution. Hence, the peak 
stress at the origin also indicates that the conical initial velocity distribution causes 
larger responses of the plate than either the spherical or cylindrical distributions. Since 
it is known that classical plate theory predicts a spurious response at all points of the 
plate immediately after application of a sharp impulse (ref. 3), the irregular stress var- 
iations occurring between t equal to 0 and 0.1 are thought to be mathematical anomalies 
rather than physical phenomena. 

The results for the displacement and stress at the origin for the three equivalent 
velocity distributions indicate, as one would expect, that a plate resists a more distributed 
impulse load better than a concentrated one. 


, 1/2 


£(2l. 

v cyl 


In all three cases a distinct peak occurs and 


APPLICATIONS 


General Momentum Impulse 

Frequently in design, the minimum plate thickness to prevent failure due to an 
impulse is required. Such a design equation can be obtained from the preceding analysis 
for a general, axisymmetric, monotonically decreasing momentum impulse on an infinite 
plate. 

If Mv'(? 7 ) is the momentum per unit area imparted to the central portion of rad- 
ius a of an infinite plate where M is the magnitude at 77 = 0 and v’(r)) is a dimen- 
sionless function specifying the distribution, the initial velocity of the plate is 


v(v) = 

v(tj) = 0 


(27) 


(0 < 77 < 1 ) 

(1 < 7J < °°) 
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Substitution of j^-v' ( 77 ) for 


v(tj) 


in equation (13) gives 


u(t) = 



( 28 ) 


The M/ph factors are specific values of V cy j, V co , and V sph of the previous 
example velocity distributions for cylindrical, conical, and spherical momentum distri- 
butions, respectively. The thickness of the plate which will just resist failure is of main 
interest. If one assumes that the response of the plate is bending and elastic until failure, 
the maximum stress at failure will be some critical stress a cr and the plate failure 

thickness some hf. Let (3L ) be the maximum value of oL where 

\ H /max p 





du 


Equation (28) then becomes 


or when solved for the failure thickness, 

h = (°p)max /3(1 
1 O nr . \ 1 


1/2 



max 


+ m)E 
-M P 



(29) 


(30) 


Hypervelocity Impact 

Spacecraft designers need design equations to specify the thickness of both single- 
and double-walled space vehicles to prevent penetration by particulate debris (meteoroids) 
while in space. Equation (30) can be applied to the impact of a cometary meteoroid with 
a single-walled space structure and to the impact of the spray emanated from the bumper 
with the main wall of a double-walled structure. However, the discussion here will be 
restricted to cometary meteoroid impact on single-walled space structures since the 
momentum per unit area distribution of the meteoroid -bumper spray cannot be defined at 
this time in terms of bumper and meteoroid parameters. 

The assumptions made about the phenomena are: momentum exchange is the pri- 
mary mechanism of interaction; the time of exchange is instantaneous; and the momentum 
of the impacting mass is negligible with respect to the momentum of the wall after impact. 
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For axisymmetric cometary meteoroids, 


M = p pip V p (31) 

where p , l , and V are the meteoroid density, length, and velocity, respectively. 

Ir ir Ir 

For the simple axisymmetric examples considered here the meteoroid mass density p i 
can be expressed as 


P 


P 


f p m p 

7Td 2 ip 


(32) 


where fp is a geometric constant. For cylinders, cones, and spheres, fp has values 
of 4, 12, and 6, respectively. Equations (31) and (32) substituted into equation (30) yield 



f 3(l + p) E\ 
\ 1 - P P) 


1/2 mpVp 

d 2 


(33) 


where 



is a factor depending only on the shape of the meteoroids. 


Low -density projectiles are difficult to launch in the laboratory; consequently, 
there are no single -wall data available which simulate the conditions assumed in this 
analysis. However, data are reported in reference 4 where plastic film disks (ref. 5) of 
different diameters impacted aluminum plates. The failure thicknesses for both projec- 
tiles were shown to vary as projectile velocity to the first power. However, if the data 
are plotted as failure thickness hf, as a function of total projectile momentum per unit 

mpVp 

area — as in figure 4, this relationship is also approximately linear. 
d z 

Values of f cfp\ can be obtained from figure 3 by noting that the curves actually 

' ' O ^ 

plotted in the three cases are 3 ? co (t), ^^^(r), and cr C yi(T). The curves were plotted 

in this manner to give the plates the same total initial momentum. One needs only to 
estimate the peak values in the figure and multiply by 1/3 and 2/3, respectively, to arrive 
at (^co)max = °’ 264 ’ ( 5 sph) max = °- 408 > and (^yl^a* = Substitution of these 
values of (^p) max along with the appropriate values of f p into the shape factors 
yields 
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(34a) 



^co 


3.17 


(34b) 


^sph " 


2.45 

n 


(34c) 


Thus equations (33) and (34) indicate that a single -walled structure of a space vehi- 
cle must be made thicker to prevent failure due to impact by a conical meteoroid than 
by either a spherical or cylindrical meteoroid when all three have the same momentum 
and diameter at impact. 

Richard Madden, Langley Research Center, 1967, employing similar methods, 
obtained an identical expression (unpublished) for equations (33) and (34a) for the par- 
ticular problem of a right circular cylinder. 


CONCLUDING REMARKS 


Linear, elastic, small -deflection plate theory was used to study the response of an 
infinite plate to an axisymmetric initial velocity distribution over a central region. The 
stress at the origin was reduced to a single definite integral which depends upon only the 
initial velocity distribution and the deflection was reduced to a Hankel inversion integral 
for three example distributions. 

Cylindrical, conical, and spherical initial velocity distributions, which gave the 
plate the same total initial momentum, were studied to demonstrate the application of the 
equations and to compare the effects of varying the distribution. The analysis indicated 
that the conical distribution caused larger deflections and stresses in the plate than either 
the spherical or cylindrical distribution. 

A simple relationship was derived for the failure thickness of an elastic plate 
responding to an axisymmetric initial momentum, monotonically decreasing in the radial 
direction. The failure -thickness relationship was shown to correlate with available exper- 
imental hypervelocity impact data where the initial momentum was induced by disk pro- 
jectiles of various diameters. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., February 4, 1969, 

124-08-01-21-23. 
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APPENDIX 


EVALUATION OF INTEGRALS 

Three integrals which appear in the text are evaluated in this appendix. The three 
integrals are evaluated by considering them to be Mellin transforms of particular func- 
tions and by selecting the appropriate value for the transform parameter. 

By definition the Mellin transform of a function is 


The integral 


p °° 

g(s) = \ f(u)u s_ l du 
J 0 

j* V 2 J 0 (PV) d?7 


which appears in equation (19a) can be considered as the Mellin transform of 


(Al) 


Hv) = J u (pv) 

i(T]) = 0 


(0 < rj < 1) 
(1 < 77 < °°) 


(A2) 


with v equal to zero and s equal to 3. The transform of equation (A2) is given in 
reference 6 (p. 326) as 

F (s±JL ;u+ 1; _I p 2) 

2^(3+ ^)r(^+ 1) \ 2 2 4 / 


which yields 


The integral 


So dl > - J i F s(fs 1 >fs-J p2 ) 

77(1 - 77 2 ) 1 / 2 J 0 (p??)dTJ 


(A3) 


(A4) 


(A5) 


which appears in equation (23a) can also be evaluated as a Mellin transform. On 
page 327 of reference 6 the Mellin transform of 


17 



APPENDIX 


i(n) = (1 - tj2) J v (pr)) (0 < 17 < 1) 

£(77) = 0 (1 < T ] < «>) 


(A6) 


is given as 


p^b(x + i,|b + y ) 

2 W_1 r(^ + l) 


p I s + , 1 s + ^ , 1 , > . 1 „ 2 

’ y+ + 1 + 4 p 


(A7) 


provided Re A > -1 and Re s > -Re v. With v = 0, s = 2, and A = , equation (A5) 

is evaluated as 


7j(l - 77 2 ) 1/2 JqCpt?) d7 7 = | 1 F 2 (l;l,|;-ip 2 


(A8) 


The identical numerator and denominator parameters of the hypergeometric function can- 
cel and leave ^ p^j . This hypergeometric function is related to the Bessel func- 

tions by the relationship (ref. 7, p. 100) 


J >) = 



r(v + i) 


o F i 




(A9) 


By use of this relationship, the hypergeometric function can be expressed as 


and the final result as 



p- 3/2 j 3 / 2 (p) 


The integral 



(A10) 


(All) 


(A 12) 


which appears in equation (25a) requires a simple change of variables to convert it to a 

Mellin transform. Let x = rp and b = — ; then the integral becomes 

4 r 
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— C (1 - x)*/ 2 cos(bx) dx 
8r Jo 

On page 320 of reference 6, the Mellin transform of 


(A 13) 


f(x) = (1 - xj^cosfbx) (0 < x < 1; Re v > 0) 


(A 14) 


f(x) = 0 


(1 < x < °°; Re v > 0) 


is given as 


| B(s,^)[^F 1 (s;s + iqib) + 1 F 1 (s;s + iq-ib)J 


(A15) 


for Re s > 0. By letting v be 3/2 and s be 1, the desired integral becomes 


§* (1 - x) 1/2 cos(bx) dx = i |^ 1 F 1 (l;|;ib) + ^ (l;|;-ibj (A16) 


The confluent hypergeometric function is the power series 

t2 


,F Aa-B-Q = 1 + SJL + °l a + 1 lll + . . 

! ^ /3 1 '. / 3(/3 + 1 ) 2 ! 

Hi 

(/3) n n! 


V ^n £ n 

1 F 1 (of;/3;^) = £ — i ~ 


(A17) 


n=0 


where Pochhammer's notation used for the numerator and denominator parameters is 


(a:) n = a(a + l)(a! + 2) . . . (a + n - 1) 


(A 18) 


and 


H = 1 


Adding the confluent hyper geometric functions of equation (A 16) term by term eliminates 
all the imaginary terms. In general, then 


oo 

+ jF^a^j-i!) = 2 ^ (-l) n 


n-0 


Hn^ 

^>2n (2n): 


(A19) 
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These results give 


1_ 

4t 




(l) 2n (4r)- 2n 

(I) <2n): 

XZ, 2n 


(A20) 


which reduces to 


h £ 1,(1 - * )1/2cos {£h “ -4- 1 

n=0 


(-l) n (4r)-^ 2n+1 ^ 

r( 2 „ ♦ |) 


(A21) 


An asymptotic expansion of the series for very small r is easily obtainable by adding 
the asymptotic expansions of the two confluent hypergeometric functions in equation (A 19). 
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(a) t = 0.15. 

Figure 2.- Comparison of the plate's deflection at various nondimensional times for the three initial velocity distributions 




Figure 2.-- Continued. 





Figure 2.- Concluded. 




Figure 3.- Variation of the maximum stress at the origin with time for the three example initial velocity distributions, 
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Figure 4.- Experimental variation of plate failure thickness with projectile momentum per unit area. Targets, 2024-T3 aluminum. 

Data obtained from reference 4. 
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